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Abstract 

We consider two-dimensional (2D) localized vortical modes in the three-wave 
system with the quadratic (y^) nonlinearity, alias nondegenerate second-harmonic¬ 
generating system, guided by the isotropic harmonic-oscillator (HO) (alias parabolic) 
confining potential. In addition to the straightforward realization in optics, the 
system models mixed atomic-molecular Bose-Einstein condensates (BECs). The 
main issue is stability of the vortex modes, which is investigated through com¬ 
putation of instability growth rates for eigenmodes of small perturbations, and 
by means of direct simulations. The threshold of parametric instability for 
single-color beams, represented solely by the second harmonic (SH) with zero 
vorticity, is found in an analytical form with the help of the variational approxi¬ 
mation (VA). Trapped states with vorticities (+1, —1, 0) in the two fundamental- 
frequency (FF) components and the SH one [the so-called hidden-vorticity (HV) 
modes] are completely unstable. Also unstable are semi-vortices (SVs), with 
component vorticities (1,0,1). However, full vortices, with charges (1,1,2), 
have a well-defined stability region. Unstable full vortices feature regions of ro¬ 
bust dynamical behavior, where they periodically split and recombine, keeping 
their vortical content. 

Keywords: vortex, second-harmonic-generation, parametric instability, azimutal 

instability, quadratic nonlinearity 
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1 Introduction 

The fundamental significance of the quadratic (x <:2> ) nonlinearity in optics, in- 
eluding its use for the creation of solitons, is well known IB-01. 0- In particular, 
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the nonlinearity opens the way to the making of two- and three-dimensional 
(2D and 3D) solitons, because, unlike the Kerr (cubic) terms, the quadratic ones, 
which couple the fundamental-frequency (FF) and second-harmonic (SH) com¬ 
ponents of the optical fields, do not give rise to the collapse in the two- and 
three-dimensional (2D and 3D) space [§,, which is a severe problem for multi¬ 
dimensional solitons in Kerr media [4 . Stable (2+l)D beams propagating in 
quadratically nonlinear optical media have been demonstrated in the experi¬ 
ment 0, and stable spatiotemporal “light bullets” have been predicted in these 
media too [TO] . In the experiment, fully self-trapped “bullets” have not been 
reported yet, the closest result being a spatiotemporal soliton self-trapped in 
the longitudinal and one transverse directions, while the confinement in the 
additional transverse direction was imposed by the guiding structure mmi]. 

A challenging issue for the (2+l)D setting is the search for conditions secur¬ 
ing the stability of vortical solitons. In contrast to their fundamental counter¬ 
parts, vortex solitons supported by the quadratic instability in the free space 
are always unstable against azimuthal perturbations, which tend to split the 
vortex into a set of separating fragments. For the degenerate quadratic system, 
with a single FF component, the splitting instability was predicted theoretically 
m-m and demonstrated in the experiment p~8] . The same instability occurs 
in the framework of the full three-wave system, which includes two distinct FF 
components, which represent mutually orthogonal polarizations of light [19, 20]. 

Vortex solitons are also known as solutions of the 2D nonlinear Schrodinger 
(NLS) equation with the self-focusing cubic term [2T . They too are subject to 
the azimuthal instability, which is actually stronger than the collapse instability 
driven by the cubic nonlinearity [4]. The 2D NLS equation models not only 
light beams in bulk media with the Kerr nonlinearity, but also the mean-field 
dynamics of atomic Bose-Einstein condensates (BECs), trapped in the form of 
“pancakes” by the confining potential [22]; in the latter case, the NLS equa¬ 
tion is usually called the Gross-Pitaevskii (GP) equation. In the context of 
the NLS/GP equation, a practically relevant solution of the instability problem 
was elaborated: both fundamental and vortical solitons, with topological charge 
m = 0 and 1, respectively, can be stabilized by 2D harmonic-oscillator (HO) 
confining potentials. It is well known 1231 - 122 ] that the HO potential renders the 
fundamental solitons stable, against small perturbations, in the entire region of 
their existence. The vortices with m = 1, trapped by the same potential, are 
stabilized in ~ 33% of their existence region (in terms of their norm). Addition¬ 
ally, an adjacent region of width 10% supports a robust dynamical regime 
featuring periodic splitting and recombination of the vortices, which keep their 
topological charge m As concerns the BEC realization, stabilization of self- 
trapped semi-vortex modes in the free 2D space was recently demonstrated in a 
two-component system with the Kerr nonlinearity and linear spin-orbit coupling 
between the components [30 . 

The effective 2D trapping potential can be also realized in optical waveg¬ 
uides, in the form of the respective profile of the transverse modulation of 
the local refractive index [5]. This circumstance suggests a natural possibil¬ 
ity for the stabilization of (2+l)D vortex solitons in the x^ medium by means 
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of the radial HO potential. For the degenerate (two-wave) x^ system, alias 
the Type-I second-harmonic-generating interaction IB-0, this possibility was 
demonstrated in Ref. m, and then extended for 3D spatiotemporal vortex 
solitons in Ref. [32] . The subject of the present work is to investigate this stabi¬ 
lization mechanism for the nondegenerate three-wave system, with the Type-II 
quadratic interaction, where the results demonstrate essentially new features in 
comparison with the degenerate model. 

A feasible approach to the making of the optical medium combining a nearly- 
parabolic profile of the refractive index and x® nonlinearity is the use of a 2D 
photonic crystal, which can be readily designed to emulate the required index 
profile, while the quadratic nonlinearity may be provided by poled material 
filling voids of the photonic-crystal matrix 0, 0. In fact, the effective ra¬ 
dial potential provides for sufficiently strong localization of the trapped modes, 
therefore the exact parabolic shape of the radial profile is not crucially important 

m- 

Essentially the same system of GP equations for atomic and molecular wave 
functions models the BEC in atomic-molecular mixtures [33]-[36]. In that case, 
two different components of the atomic wave function pertain to two different 
atomic states of the same species, while the quadratic nonlinearity accounts 
for reactions of the merger of two atoms into a molecule, and splitting of the 
molecules due to collisions with free atoms. Accordingly, the predicted mecha¬ 
nism of the stabilization of three-component vortex solitons trapped in the HO 
potential can also be realized in the BEC mixture. 

The paper is organized as follows. The model, based on the system of three 
propagation equation coupled by the x^ terms, is introduced in Section II. 
In Section III, we introduce four types of three-wave modes considered in this 
work: the SH-only (single-color) one, with zero FF components; modes with 
the hidden vorticity (HV, which means opposite vorticities, S = ±1, in the two 
FF components, and 5 = 0 in the SH component); semi-vortices, with S = 0 
in one FF component and S' = 1 in the other one and in the SH wave (cf. 
the above-mentioned semi-vortices in the two-component system with the cubic 
nonlinearity and linear spin-orbit coupling between the components [30]); and, 
finally, full vortices , with S = 1 in both FF components, and S = 2 in the SH. 
In Section IV we report simple but nontrivial analytical results, which predict 
the threshold of the parametric instability for the fundamental and vortical 
single-color states by means of the variational approximation (VA). In Section V 
we formulate the eigenvalue problem for small perturbations around stationary 
solutions, which determines their stability. The numerical results for modes of 
all the types are reported in Section VI. First, we identify the stability threshold 
for the fundamental and vortical single-color trapped modes, and compare these 
results with the above-mentioned analytical predictions provided by the VA. 
Next, the results are summarized for the HV and SV modes, which are found 
to be unstable. Finally, the stability region is presented for the most general 
three-wave states, with vorticities 1 in both FF components and 2 in the SH. 
In addition, unstable vortices may exist in a quasi-stable dynamical regime, in 
the form of periodic splitting and recombination of the semi-vortex. The paper 
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is concluded by Section VII. 


2 The model 

The copropagation of two components of the FF wave, with amplitudes u {x, y , z) 
and v (x, y , z), and the SH wave, with amplitude w (x, z), in the bulk x ^ waveg¬ 
uide obeys the system of three coupled equations written here in the normalized 
form [213]: 


iu z = — - V 2 u + v*w + -^t 2 r 2 u + Qu, (1) 

iv z = — + u*w + ^Q 2 r 2 v — Qv, (2) 

2 iw z = — — qw + uv + 2Q, 2 r 2 w, (3) 


where V 2 = d 2 /dx 2 -\-d 2 jdy 2 is the transverse diffraction operator, the asterisk 
stands for the complex conjugate, Q is the birefringence coefficient, q is the 
FF-SH phase mismatch, and Q 2 the strength of the HO trapping potential. 
In the most general case, stationary vortex solutions of this system, with two 
independent FF propagation constants k u and k v , can be looked for as 

{u,v} = {U 0 (r)e iS * e+ik '' z , V 0 (r)e iS '’ e+ikvZ } , (4) 

w = Wb(r)e < ( 5 “ +5 ^ +< ( fc “ +fc «^, (5) 


where S u and S v are integer vorticities of the two FF components, and Uo(r), 
Vo(r), and Wo(r) are real radial functions with asymptotic forms 

{Uo(r), Vo(r), Wo(r)} ~ r { |-S' u | ,|5^|,|5 u +5^|} at r ^ 0 , 

{Uo(r), Vo(r)} ~ exp 5 Wo(r) ~ exp (— fir 2 ) at r oo. (6) 

The exact equality of the propagation constant and vorticity of the SH compo¬ 
nent in Eqs. Q and (J5| to sums of the same constants of the FF components is 
imposed by the coherent coupling between the FF and SH components in Eqs. 
0, (J2|, and 0. 

In this work, the analysis is carried out, chiefly, for Q = 0, in which case 
the symmetry between the two FF components suggests to consider three-wave 
states with k u = k v . For Q ^ 0, rescaling of Eqs. 0-0 makes it possible to 
fix Q = 1. Families of three-wave vortices with Q = 1 are considered too, at the 
end of the paper. 

Equations 0 -© can be derived from the Lagrangian, 

L = ( f T [(uu* z - u*u z ) + i ( vv* z - v*v z ) + 2 {ww* z - w*w z )\ 
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dxdy — H, 

(7) 


where the Hamiltonian is 


H = 11 \ (|Vu| 2 + |V«| 2 + | Vw| 2 ) + Ul 2 r 2 (|u| 2 + \v\ 2 + 4M 2 ) 


+Q (\u\ 2 — \v \ 2 ) — q\w\ 2 + ( uvw * + u*v*w)] dxdy. (8) 


( 8 ) 


The Hamiltonian is the dynamical invariant of Eqs. 0-®> along with the 
total power (field norm), 



and the total angular momentum, 


M = i [u* ( yu x — xu y ) + v * (yv x — xv y ) + 2 w* ( yw x — xw y )\ dxdy , (10) 



which is a real integral quantity, even if it formally seems complex, (it can be 
checked by means of the integration by parts that the M* = M). 

3 Stationary solutions 

3.1 Single-color (second-harmonic) states 

The simplest stationary state in the present model is the one with unexcited FF 
components, i.e., Uo = Vo = 0 in Eq. <0L and Wo(r) from Eq. © obeying the 
linear equation, which follows from Eq. 



( 11 ) 


where S w replaces (S u + S v ), see Eq. ©• This equation is tantamount to the 
radial Schrodinger equation for 2D quantum-mechanical states with azimuthal 
quantum number S in the HO potential. The respective solutions to Eq. © 
are 


W 0 = Brl s » l exp(-Or 2 ), 
k = —12 (1 + ISwl) /2 + q/4, 


( 12 ) 

(13) 


where B is an arbitrary constant. The norm © of this state is 



( 14 ) 
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3.2 Hidden-vorticity (HV) states 

Solutions for 2D HV modes are defined by the vorticity set (S u , S V ,S W = S u + S v ) 
(+1, —1,0) of the three components in Eqs. © and ©• Accordingly, the solu¬ 
tions are looked for as 

{u, v} = {Uo(r)e iSe , V 0 (r)e~ iSe } e ikz , (15) 

w = W 0 (r)e 2ikz , (16) 

where real functions Vo, Vo, and Wo satisfy the following radial equations (where 
the birefringence terms are kept, for the time being): 

- kU 0 = -^^ + l -U' 0 -^U^+V 0 Wo + \ti 1 r 2 U 0 + QUo, (17) 

~kV 0 = ~l(v£ + ^-^V 0 ^+UoW 0 + ^ 2 r 2 Vo-QV 0 , (18) 

—4fcWo = ~\WS + ^w£j-qWo + UoVo + 2to 2 r 2 Wo. (19) 

Previously, the HV concept was realized for the three-wave x ^ system in 
the free space El, as well as for various systems of coupled continuous m-m 
and discrete [42j 03] NLS/GP equations with the cubic nonlinearity, including 
the system supported by the HO trapping potential [41 . In work [44], it was 
demonstrated that the HV essentially suppresses the modulational instability 
of 2D ring solitons in a two-component system with saturable nonlinearity, in 
comparison with their counterparts carrying explicit vorticity. However, trapped 
vortical modes with HV were not previously studied in three-wave systems. 


3.3 Semi-vortices 


The other type of the vortical mode, with zero vorticity in one FF component 
(hence they are called semi-vortices, as said above), is defined by the set of 
(S u , S v , S w ) = (1,0,1), i.e., solution ansatz ©,© reduces to 

{«,«} = {U 0 (r)e ie+ik ' z ,V 0 (r)e ik2Z }, (20) 

w = W 0 (r)e ie+i(kl+k2)z , (21) 


where the real radial amplitudes obey the following equations, cf. Eqs. 

(l8|, ©: 


— kiUo 


-k 2 V 0 


-2 (h + k 2 ) Wo 


(vs + Vo - >o) + V 0 w 0 + ^n 2 r 2 U 0 + QU 0 , 
-\ [Vo + Vo) + U 0 Wo + VW 0 - QV 0 , (22) 
(w" + )w ( ( - 4 Wo ) _ qWo + UqVq + 2 n 2 r 2 W 0 . 
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In the limit of small amplitudes, solutions of the linearized version of Eqs. (22) 
can be found in the form of 


U 0 = U^rexp (—Qr 2 / 2) , 

V 0 = W 0 (0) exp (-Or 2 /2) , (23) 

Wo = Wg 0 Vexp (—f2r 2 ) , 

where the propagation constant must satisfy, respectively, the following relations 
for the components U,V,W: 


ki = —2fi - Q, 

k 2 = —Q + Q, (24) 

k\ H - k 2 = + qj 2, 


From Eqs. (24) it follows that, in the small-amplitude limit, the semi-vortex 
mode exists at the single value of the mismatch, q = —2f2. Below, it is shown 
that solutions for nonlinear semi-vortex modes lift this constraint. 


3.4 Full vortices 


In fact, the HV modes and semi-vortices are shown below to be completely unsta¬ 
ble (modes of the latter type may be replaced by robust dynamical states which 
keep the same vortical structure). A stationary topological state which features 
a well-defined stability region is the full vortex , with the set of (S u , S v , S w ) = 
(1,1,2) in Eqs.j4]) and ([5|. The respective system of stationary equations is 
[cf. Eqs. @-([I^and @]: 


-hUo 


-k 2 V 0 


-2 (A* + k 2 ) Wo 


-4 (u{{ + To - » + VoWo + \n 2 r 2 Uo + QUo, 
-\ (V 0 " + T 0 ' - + U 0 Wo + ifi 2 r 2 F 0 - QV& 5) 

\WS + -W^ - 4 W 0 ) - qWo + UoVo + 2 fi 2 r 2 W 0 . 


4 The variational approximation (VA) for the 
onset of the parametric instability of the single¬ 
color modes 

The VA is a natural approach for predicting the shape of nonlinear modes, 
including vortical ones supported by the nonlinearity mm- In the present 
context, the VA can be used for the prediction of the instability threshold of 
the trapped single-color (SH-only) states (the VA may be developed for other 
situations too, but it then turns out to be rather cumbersome). 
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For this purpose, it is necessary to consider a perturbed version of the single¬ 
color stationary state, which includes the FF components too. In particular, in 
the case of zero vorticity of the SH component and vorticities ±S of the FF 
perturbations, the Lagrangian of the respective stationary equations ([T7|)-([T9|) 
can be reduced to its radial part: 



+k (U 2 + V 2 + 4W 2 ) + I n 2 r 2 (U 2 + V 2 + 4W 2 ) + Q ( U, 2 - V 2 ) - qW 2 

+2U 0 V 0 W 0 }. (26) 


To develop the VA, the simplest Gaussian ansatz may be adopted, assuming 
equal amplitudes ( A ) and widths of fields Uo and Vb: 

U 0 (r) = V 0 (r) = Ar s exp (-Or 2 /2) , W 0 (r) = B exp (-Or 2 ) (27) 


(S > 0 is defined here). The symmetry between the FF components implies that 
we should set Q = 0, otherwise the symmetry will be broken by the birefringence. 
Setting Q = 0 limits the consideration to a particular case of the generic three- 
wave system, but even this particular case is a nontrivial one, as it has no 
counterpart in the usually considered degenerate two-wave system, in which the 
VA was applied before to the description of trapped vortices [31] . The widths of 
all the components in ansatz (27) are not treated as variational parameters, but 
are rather taken as per wave functions of states trapped in the 2D harmonic- 
oscillator potential, cf. Eq. (12). 

The objective of the use of the VA in the present context is to predict a point 
at which a solution with an infinitely small A appears, against the background 
of the fundamental single-color (SH-only) state, given by Eqs. (12), (13) with 
S w = 0 and arbitrary amplitude B. The appearance of this mode signals the 
onset of the parametric instability of the latter state [3F. The substitution of 
ansatz (27) into radial Lagrangian (26) yields the following reduced Lagrangian, 
in which we drop terms that produce no contribution in the subsequent analysis: 

(28) 


^rad 


n s 


1 + s + Ii +2 ~ {1+S) fi 


This Lagrangian gives rise to the Euler-Lagrange equation, dL^/d (A 2 ) = 0. 
The parametric-instability threshold, corresponding to the emergence of the 
solution with an infinitesimal A 2 > 0, corresponds to setting A 2 = 0 in the 
resulting equation: 


B = -2 1+s [(1 + S) n + k] = —2 1+s 


- + s + 4 


(29) 


where expression (13) was substituted for k. The respective critical total power 
(norm) of the trapped single-color state is given by Eq. (14), 


N c = 2 1+2S ^ [(l + 2S)fi+| 


2 


(30) 



















Thus, the fundamental single-color state is stable at N < N c , and unstable at 
N > N c . 


Note that, as it follows from Eq. (30), the parametric instability of the zero- 


vorticity single-color state is dominated by the vortical FF perturbations with 
S = 1 (i.e., S = 1 give rise to N c lower than its counterpart produced by the 
zero-vorticity perturbations, with S = 0) in an interval of negative values of the 
mismatch, (14/3) Q < — q < lOfh 

It is also relevant to develop a similar analysis for the onset of the parametric 
instability of the vortical single-color trapped state with even vorticity S w > 0, 
for which Wo(r) = Br Sw exp (—fir 2 ), cf. Eq. (12), and S is replaced by S w /2 in 
the FF components of ansatz (|27|). In this case, the reduced radial Lagrangian 


(28) is replaced by 


bfd” } = 0W2)! 


Afi + ^ + *bs » A2J? 

n s w /2 2 ft; (2ft) S ” +1 ‘ 


(31) 


Then, the Euler-Lagrange equation, /d (A 2 ) = 0, yields the value of 

amplitude B , the substitution of which into expression ( [l4| ) gives the instability 
threshold for the single-color vortex mode: 


N (S W ) = 2 1+S W [(S w /m 2 ^ ( Q , Q y 

c S w \ fiv + 2/ 


(32) 


to numerical findings below, see Fig. |TJ 


The analytical predictions (30) and (|32j) produced by the VA are compared 


5 Linearized equations for small perturbations 

The stability of the modes under the consideration is the central issue of the 
present work. It is addressed via computation of eigenvalues for small pertur¬ 
bations. In the general case, perturbed solutions are introduced as 


u(z,t) 

= e ikz 

'u 0 {r) + £e Xz Ul {r) + ee x ' z ul{r) 



v 0 , r) 

= e ikz 

t>o(r) +ee Xz vi{r) + ee x z vl(r) 


(33) 

w (z, r) 

_ ^2 ikz 

wo(r) + ee Xz wi(r) + ee x z w{( r) , 



where k is the propagation constant of the stationary solution (assuming here 
that k is the same for that both FF components), which is represented by 
(complex) functions ^o( r ) 5 ^o( r )> wo(r), infinitely small 5 is an amplitude of 
the perturbation, A is the instability growth rate sought for, which is, generally, 
complex too [the instability takes place if there exists, as usual, at least single A 
with Re(A) > 0], the asterisk stands for the complex conjugate, and u i(r), fti(r), 
^i( r ) 5 fh( r ) 5 ^i( r ) are components of the perturbation eigenmode, which 
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obey the system of linearized equations: 


(iX — k)ui = -^V 2 rq + v^wi + + ifl 2 r 2 rq + Qui, 

(~i\ — k)ui = — ^V 2 tq + + ^Cl 2 r 2 ui + Qiii, 

(iA — fc) t?i = —^V 2 fi + u^wi + W 0 U 1 + if^ 2 r 2 fi — Q^i, 

(—iA — fc) vi = —^V 2 hi + ^ 0^1 + Wo u i + ^fl 2 r 2 vi — Qvi, 

(2iX + q — 4 k) w\ = — ^S7 2 w\ + uqV\ + v^ui + 2Q 2 r 2 wi, 

(—2iX + q — 4 k) w\ = — ^V 2 uq + UqVi + VqUi + 2Q 2 r 2 wi. (34) 


For the unperturbed solution taken as per Eqs. (15) and (16), it is easy 
to see that Eqs. (34) admit self-consistent perturbation modes in the following 
form: 


iq(r, 0) = £/i(r)e i( -^ +p)0 , fq(r, 0) *= tfi(r)e i(5 " +p) *, 

vi(r,0) = Vi(r)e^-^ + ^, ®i(r, 0) = Vi(r)e^+^, 

wi(r,0) = ^(r, 0) = W 1 (r)e i ( s * +s '’ +p)e , (35) 

where p is an arbitrary integer (azimuthal index of the perturbation). The 
substitution of expressions ( [35] ) into Eq. p4| ) leads to the following eigenvalue 
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system which should determine A for each integer p : 


= + ( S ; 2 +y)2 ) tr, + H(W + IVoMf-, 

+^VC/i + QC/i, 

(-a -1) ft = -i (^ + ~ - ) ft + HOW + »'„(r)n 

+^nWi + Qt7 lt 

(iA - t) Vi = -i + H - ( ~ S lt r}2 ) K + ft(W + IVoMft 

+ hr 2 r 2 l/ 1 - QV U 

( _a - i) a = -1 (Ji + I ± - 2^1!) a + v„(rm+w^m 


+ -n 2 r 2 V 1 -QVu 

. A . _ . Txr 1 f d 2 Id (~S U — S v + p ) 2 \ _ _ _ 

(2,A + <, - 4fc) in, = -- ^ - / j >ft + ft W 

+Vo(i')!7i+2nVW' 1 , 

+U 0 (r)V 1 +V 0 (r)U 1 +2n 2 r 2 W 1 . (36) 


Solutions U \, U\ , Vi , Vi of Eqs. (|36j) must exponentially decay at r —>• oo, and 
behave as r \p ±s u,v\ at r —>> 0, and Wi, IVi must exponentially decay too at 
r —)► oo, and go as r \p±(Su+s v )\ a t r —)> 0, cf. Eq. |6|. 


6 The stability of the single-color modes 


First, we address the stability of the SH-only vortex, which is given by Eqs. 


(12) and (13), against small perturbations seeded in the FF components, which 


are taken as 


u(z, r) = ee ikz+Xz Ul (r)e ipe , 

v(z, r) = ee <fc*+A** u .( r ) e i(s»-p)« > ( 37 ) 


according to Eqs. (33), ( [35] ) . The perturbation eigenfunctions, 'Ui(r) and v>\ (r), 
satisfy the following linearized equations, which are a particular case of Eq. 
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( 36 ): 


1 ( d 2 Id P 2 \ 

+£?r‘ s “’e _nr + if2 2 r 2 wi, 

1 / d 2 Id -P) 2 ' 

(_ tA _ A + Q )t;i - 


Vi 


+Br Sw e~ nrQ 


u\ + i^ 2 r 2 t’i. 


The objective is to find, for given S^, a critical (minimum) value of amplitude 
B in solution (Jl2|, and, accordingly, the minimum value of the integral power 
(14), at which Eq. © starts to produce eigenvalues with Re(A) ^ 0, i.e., the 
single-color mode becomes unstable against the FF perturbations. Then, we aim 
to explore the evolution of unstable modes by means of direct simulations. To 
address these problems, we employed numerical techniques and grids similar to 
those used in Ref. m for finding the eigenvalues and running direct simulations. 

The results for trapped single-color (SH-only) modes are collected in the left 
panel of Fig. [l] while the right panel illustrates the evolution of an unstable 
mode with S w =0, in terms of the power exchange between the three compo¬ 
nents (the total power is conserved, as it should be). Here and below, numerical 
results are displayed for the zero birefringence (Q = 0, unless it is specified 
otherwise) and trapping-potential strength ft = 1. 

The left panel of Fig. [l] demonstrates that the VA-predicted instability 
thresholds ( |3()| ) and ( 30 ) approximate their numerically found counterparts well 
enough. In particular, the VA exactly predicts that (for ft = 1 and Q = 0), the 
instability threshold vanishes (N c = 0) at qm —2. 

It is interesting to consider the (in)stability of the single-color trapped vor¬ 
tices with S w = 1 and other odd values of the vorticity, as the respective per¬ 
turbations in the two FF components cannot be arranged symmetrically, unlike 
those considered above for even values of S w . The threshold value of the norm 
for S w = 1 is shown in the left panel of Fig. [lj and the development of the 
respective instability is illustrated by Figs. © and ©. it can be concluded 
that the instability does not destroy the vorticity (topological charge) of the 
SH component, while the SH and FF waves exchange their angular momenta, 
keeping virtually equal integral powers. 


7 Numerical results for three-wave vortices 
7.1 Hidden-vorticity (HV) modes with S UfVtW = (+1,—1,0) 

Families of HV states, defined as per Eqs. <© and @, are characterized 
by dependences between the propagation constant and integral powers of the 
two FF components and the SH one, see Eq. ©• Typical examples of such 
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z 



Figure 1: Left: Numerically found critical powers, N c , of the single-color 
(SH-only) states, with vorticities S w = 0, 1 , 2, beyond which the states become 
unstable against the excitation of the FF. For S w = 0 and 2, the analytical 
results, produced by the variational approximation [see Eqs. (30) with S = 0 
and (32) with S w = 2], are displayed too. Right: The evolution of integral 
powers of fields u , v and w for a perturbed unstable state with S'™ = 0, g = 0, 
and total power N = 7.6. In both panels here, and in figures following below, 
numerical results are presented for Q = 1 and Q = 0 (zero birefringence between 
the FF components, except for Fig. 15). 
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Figure 2: (Color online) Development of the instability of the single-color (SH) 
mode (12), for S w = 1, Q = 0, = 1, B = 2, TV = 47r. Initial random 

perturbations in the FF components, u and u, are loaded with an amplitude 
< 10 -4 . The top, middle, and bottom panels display the evolution of densities 
of the u , v (FF) and w (SH) fields, respectively. 


dependences, obtained from a numerical solution of Eqs. (pT|)-([l9|, are shown 
in the left panels of Figs. [4]and[5j respectively, for positive (q > 0) and negative 
(q < 0) values of the phase mismatch. 

The computation of the stability eigenvalues by means of Eq. (36) demon¬ 
strates, in the right panels of Figs. [4] and [5} that all the HV modes are subject 
to instability against small perturbations with various values of azimuthal index 
p. Direct simulations reveal two different scenarios of the instability develop¬ 
ment, which are displayed in Figs. [6] and [7] In the former case, the central core 
is spontaneously expelled from the vortex ring, which tends to transform itself 
into a fundamental mode. In the latter case, the vortex ring splits into two 
fragments, which is a typical outcome of the instability development of vortices 
in models with the cubic nonlinearity [25]- [28 ] . 
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Figure 3: Overall characteristics of the instability displayed in Fig. [2j Left: The 
evolution of the integral powers of the FF components, u and v. Right: The 


integral angular momentum of each component, see Eq. (10). 


7.2 Semi-vortices with S UiViW = (1,0,1) 


For families of semi-vortex solutions, a numerical solution of Eq. (22) produces 


dependences between the propagation constant and integral powers of the three 
components, a typical example of which is displayed in Fig. [8] It is concluded 
that these families are completely unstable too. 

Simulations of the instability development of the semi-vortices reveal two 
basic outcomes which initially look similar to those shown above for the HV 
modes, cf. Figs. [6]and[7| i.e., expulsion of the central core from the vortex ring, 
as shown in Fig. [9| or splitting of the vortex into a set of two fragments, see 

Fig. m 


7.3 Full three-wave vortices 

The most general three-wave vortex state is built with S UiV ^ w = (1,1,2). Fam¬ 
ilies of these solutions were found for the symmetric system, with Q = 0 and 
k u = k v = k [see Eqs. Q and ([5|] from a numerical solution of Eq. ( [25] ) . Their 
stability eigenvalues were then computed using Eq. (36), and the predicted 
stability or instability was verified by direct simulations of Eqs. 0 -®>- 

A crucial difference of the three-wave vortices from the single-color (SH- 
only), HV and SV states, which were considered above, is that the full vortices 
have a well-defined stability area. A typical vortex family and its stability are 
presented in Fig. [TT] This figure explicitly displays both the bifurcation, which 
generates the full vortex from the corresponding single-color state, with S w = 2 
(as shown above, the bifurcation simultaneously destabilizes the single-color 
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Figure 4: (Color online) Left: The propagation constant vs. the total power 
for the hidden-vorticity family with vorticities S UiViW = (+1,—1,0), D = 1, 
Q = q = 0. Right: Stability eigenvalues produced by linearized equations (36). 
The corresponding values of the perturbation azimuthal index p are attached to 
the curves. 


state), and the point of the destabilization of the three-wave vortices. 

The results of the stability analysis are summarized, in the plane of the 
phase-mismatch parameter, g, and total power, TV, in Fig. [T2| The shape 
of the stability diagram is qualitatively similar to the one which was recently 
reported, for the degenerate two-wave x ^ system, in Ref. |3l] (for the sake of 
the comparison, note that q was defined with the opposite sign in Ref. EB). 
The bottom-right area in the diagram is populated by the single-color (SH-only) 
vortices with S w = 2, the bifurcation destabilizing the single-color vortex and 
replacing it by the three-wave one occurring along the upper boundary of this 
area. 

In the “splitting” area labeled in Fig. [l2j the three-wave vortices are subject 
to an instability which splits them into a set of fragments, the number of which 
is equal to the dominant (or single) azimuthal index, p, of unstable perturbation 
modes, which is indicated in the right panel of the figure. Further, in the region 
labeled “p = 2 (only)”, the stable static three-wave vortices are replaced by a 
robust dynamical regime, in the form of periodic splittings of the vortex into 
two segments and their recombinations, as shown in Fig. 
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In the course 

of this periodic evolution, the vortical structure of the mode is conserved. A 
similar dynamical regime, in the form of periodic splittings and recombinations, 
is known in the 2D GP equation with the cubic self-attractive nonlinearity and 
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Figure 5: (Color online) The same as in Fig. [4j but for q = —4. 


HO trapping potential m- On the other hand, the instability-induced splitting 
the vortex into a set of three fragments is irreversible, as shown in Fig. [l4j 
Finally, the effect of the birefringence of the stability of the three-wave vor¬ 
tices, represented by Q = 1 in Eqs. 0-0, was briefly considered too. As 
shown in Fig. [l5j the birefringence terms make the stability area somewhat 
larger. This result can be explained by the fact that the birefringence renders 
the system less coherent, while the splitting instability of the vortices is a result 
of highly coherent x ^ interactions. 

8 Conclusion 

This work aimed to explore the possibility of the stabilization of various 2D 
three-wave modes supported by the x ^ interactions in the combination with 
the isotropic HO (harmonic-oscillator) trapping potential. The existence and 
stability of the modes is determined by powers and vorticities of the three com¬ 
ponents and the mismatch of the x ^ system (q). First, using both numerical 
computations and the VA (variational approximation), stability boundaries were 
identified for the fundamental (zero-vorticity) and vortical single-color states, 
in which only the SH (second-harmonic) component is present. On the con¬ 
trary to the usual assumption that the single-color SH modes are subject to the 
parametric instability against perturbations in the FF (fundamental-frequency) 
fields, we have found that they are stable below the respective critical values 
of the total power. Next, HV (hidden-vorticity) and SV (semi-vortex) states, 
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Figure 6: (Color online) The evolution of the two components, u , of an unstable 
hidden-vorticity mode with q = —4, norms N u , v = 25.1, N w = 52.8 (the total 
power is N = 103) and k = — 1. Density distributions are displayed for values 
of the propagation distance indicated above the respective panels. 
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Figure 7: (Color online) The same as in Fig. [6j but for q = — 4 and N UiV = 340, 
N w = 480, the total power is N = 1160, and k = 2. 


with vorticities, respectively, ±1 or 1 and 0 in the two FF components, were 
found to be always unstable. The stability region has been identified for the full 
three-wave vortices. Furthermore, adjacent to it is the region which features 
the robust dynamical regime of periodic splitting into two fragments and their 
recombination into the original vortex. 

The analysis reported in the present work can be extended in other direc¬ 
tions. One possibility is to analyze asymmetric three-wave modes, with unequal 
propagation constants in the two components of the FF field. It may also 
be interesting to construct self-trapped three-wave x ^ modes supported by a 
periodic (lattice) potential, instead of the HO, cf. Ref. [45] . Lastly, a chal¬ 
lenging possibility is to construct three-dimensional three-wave “light bullets” 
supported by the HO trapping potential, cf. Ref. [32] where this was done for 
the degenerate two-wave system. 
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Figure 8: (Color online) Left: Relations between the propagation constant and 
integral powers of the three fields for semi-vortices of the type of S u ^ VrW = 
(1,0,1), with Ll = 1 , Q = 0, q = 0. Right: The largest real parts of stability 
eigenvalues for the semi-vortices, as obtained from the numerical solution of 
Eq. ( [36] ) . Integers p are azimuthal indices of the destabilizing perturbation 
eigenmodes. 
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numerical solution of linearized equations (36). 


[30] Sakaguchi H, Li B and Malomed B A 2014 Phys. Rev. E 89 032920 

[31] Sakaguchi H and Malomed B A 2012, J. Opt. Soc. Am. B 29 2741 

[32] Sakaguchi H and B. A. Malomed 2013 Opt. Exp. 21 9813-9823 

[33] Drummond P D, Kheruntsyan K V and He H 1998 Phys. Rev. Lett. 81 
3055-3058 

[34] Heinzen D J, Wynar R, Drummond P D and Kheruntsyan K V 2000 Phys. 
Rev. Lett. 84 5029-5033 

[35] Hope J J and Olsen M K 2001 Phys. Rev. Lett. 86 3220-3223 

[36] Hornung T, Gordienko S, de Vivie-Riedle R and Verhaar 2002 Phys. Rev. 
A 66 043607 

[37] Torres J P, Soto-Crespo J M, Torner L and Petrov D V 1998 Opt. Commun. 
149 77-83 


22 








Figure 12: (Color online) Left: The stability diagram for vortices with S u , v , w = 
(1,1,2), D = 1 , Q = 0, in the plane of the mismatch parameter (g) and total 
power (TV). Right: The same, on the logarithmic scale of TV. In this panel, 
dominant values of the perturbation azimuthal index are displayed in the region 
of the splitting instability (label “only” implies that there is a single instability 
eigenmode in the respective area). 


[38] Desyatnikov A S, Mihalache D, Mazilu D, Malomed B A, Denz C and 
Lederer F 2005 Phys. Rev. E 71 026615 

[39] Leblond H, Malomed B A and Mihalache D 2005 Phys. Rev. E 71 036608 

[40] Desyatnikov A S, Mihalache D, Mazilu D, Malomed B A and Lederer F 
2007 Phys. Lett. A 364 231-234 

[41] Brtka M, Gammal A and Malomed B A 2010 Phys. Rev. A 82 053610 

[42] Kevrekidis P G and Pelinovsky D E 2006 Proc. R. Soc. A 462 2671-2694 

[43] Leykam D, Malomed B and Desyatnikov A S 2013 J. Optics 15 044016 

[44] Ye F, Wang J, Dong L, and Li Y 2004, Opt. Commun. 230 219-223 

[45] Xu Z Y, Kartashov Y V, Crasovan L C, Mihalache D and Torner L 2005 
Phys. Rev. E 71 016616 


23 






u,v z=190 u,v z=200 u,v z=210 



u,v z=230 


o 


-3 

- 3 - 2-10123 


W Z=190 


© 


-3 

- 3 - 2-10123 


W Z=200 



X 


w z=210 



x 


w z=230 


© 


-3 

- 3 - 2-10123 


Figure 13: (Color online) The evolution of the densities of the u , v, w fields in 
the regime of periodic splittings and recombinations of the three-wave vortex, 
with S u ,v,w = (1,1,2), ft = 1, Q = 0,q = 0, k = -1, N u , v = 21.4, N w = 44.4, 
the total power being N = 87.2. 
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Figure 14: (Color online) An example of the irreversible splitting of the three- 
wave vortex with S U:V ^ W = (1,1,2). Top: The evolution of densities of the FF 
fields, u and v, for d == 1, Q = 0, q =» 0, k m 10, N UjV = 1800, N w = 2142, the 
total power being N = 5741. Bottom: The same but for the SH field, w. 
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Figure 15: (Color online) The same as in the right panel of Fig. 
Q = 1 in Eqs. 


11, but with 
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